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Pressure- Temperature Phase Diagram of Golden SmS 
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We measured the thermal expansion of the valence fluctuating phase of SmS (golden SmS) 
to construct a pressure vs temperature phase diagram. The obtained phase diagram is char- 
acterized by three lines. One is a crossover line that divides the paramagnetic phase into two 
regions. The other two lines correspond to a second-order Neel transition and a first-order Neel 
transition. The crossover line appears to emerge from a tricritical point that separates the 
first-order Neel transition from the second-order one. We argue that a valence jump occurs at 
the border of antiferromagnetism. 
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1. Introduction 

The hybridization between 4/ and conduction elec- 
trons forms a heavy-fermion state accompanied by 
spin and/or charge fluctuations. Although heavy-fermion 
compounds are basically metallic, they show various 
ground states, for example, heavy Fermi liquid, magnetic 
orderings and superconductivity. Interestingly, there are 
some compounds forming an energy gap or a pseudo- 
gap at low temperatures. Typical examples include SmBg 
and high-pressure-phase SmS, which we call golden SmS 
hereafter. 1 Quantum phase transition (QPT), which is a 
phase transition at zero temperature, has also attracted 
our interest, because some exotic phenomena such as 
non-Fermi liquid and unconventional superconductivity 
appear in the vicinity of QPT. In heavy fermions, rele- 
vant energy scales are so small that QPT is easily tuned 
by rather small external parameters such as pressure and 
magnetic field. In this study, we focus on the pressure- 
induced phase transition in golden SmS. 

At ambient pressure, SmS is a nonmagnetic, ionic crys- 
tal with a small energy gap (E g ~ 0.1 eV). The total 
angular momentum arising from the Sm 2+ (4/ 6 ) config- 
uration is vanishing (J = \L — S\ = 0), where L = 3 
and S — 3 denote the orbital and spin angular momenta, 
respectively. With increasing pressure, energy gap de- 
creases monotonically and finally collapses at a certain 
pressure (P c j ~ 7 kbar at room temperature). 2,3 This 
isostructural first-order phase transition involves a va- 
lence change from divalence to mixed valence of Sm 2+ 
and Sm 3+ accompanied by a color change from black to 
golden yellow. This is the reason why the high-pressure 
phase is called golden SmS. Note that a pseudogap ap- 
pears to open in this intriguing phase. 4 When external 
pressure is further increased at low temperatures, SmS 
undergoes a second phase transition at a critical pressure 
(P c2 ~ 19 kbar at T~ 0). 

Improvements in high-pressure experimental tech- 
niques have contributed to the elucidation of the elec- 
tronic state of golden SmS: Experiments on ac specific 



heat and thermal expansion revealed a sharp anomaly 
in their temperature dependence, indicating the pres- 
ence of a bulk phase transition at pressures above P C 2- 5,6 
Nuclear forward scattering (NFS) experiments clearly 
showed that the phase transition is of magnetic origin. 7 
This means that the magnetic Sm + ionic state with the 
configuration 4/ 5 is stable at high pressures above P C 2- 
Furthermore, ac magnetic susceptibility experiments in- 
dicated that the magnetic ground state is antiferromag- 
netic. 8 

The NFS experiment showed the sudden appearance of 
an internal magnetic field when pressure increases across 
P C 2, suggesting the first-order phase transition at low 
temperatures. On the other hand, specific heat and ther- 
mal expansion experiments demonstrated the second- 
order nature of the phase transition at pressures greater 
than P C 2 . The combination of these two results suggests 
that there is a tricritical point (TCP) that separates the 
first-order Neel transition from the second-order one near 
P c2 in the pressure (P) vs temperature (T) phase dia- 
gram. However, a complete P — T phase diagram is still 
missing. In this paper, therefore, we present a P— T phase 
diagram of golden SmS using the thermal expansion tech- 
nique, in which the presence of the TCP is confirmed by 
the observation of a discontinuous volume change in the 
vicinity of P C 2 . 

From an experimental point of view, thermal expan- 
sion is a very useful tool for studying QPT, because it 
will be largely enhanced near QPT owing to a diverging 
Griineisen parameter. 9, 10 On the other hand, note that 
a pressure-transmitting medium has a crucial effect on 
physical properties including thermal expansion. In our 
previous report, pressure hysteresis was observed in the 
P — T phase diagram of SmS above 10 kbar. 6 This pres- 
sure roughly coincides with the "solidification" pressure 
of Fluorinert that was used as a pressure-transmitting 
medium, and hence the hysteresis may not be intrinsic. 
In this study, we first determined the solidification pres- 
sure of Daphne oil 7373 as a function of temperature, and 
then remeasured the thermal expansion of SmS single 
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crystals using Daphne oil to obtain the intrinsic effects 
of pressure on SmS. 

2. Experimental Results and Discussion 

2.1 Experimental methods 

SmS single crystals were grown from starting materials 
of 99.99% pure (4N) Sm chips and 6N powdered S by the 
vertical Bridgman method in a high-frequency induction 
furnace. 11 Three types of single crystals with different 
starting compositions (nominally stoichiometric, 1% ex- 
cess of Sm, and 1% excess of S samples) were prepared. 
A sample with typical dimensions of 1.5 x 1.6 x 0.45 
mm 3 was cleaved for measurements. Thermal expansion 
results were found to be only weakly sample-dependent, 
so that we report here the results of nominally stoichio- 
metric samples. 

Thermal expansion measurement under high pressure 
was carried out by the active-dummy method using 
strain gauges (SKF-5414, KYOWA). A copper block (6N 
purity, typical dimensions of 1.5 x 1.5 x 1.0 mm 3 ) was 
used as a dummy sample. We employed an ac method 
in terms of a lock-in amplifier to reduce the drift voltage 
arising from thermoelectric power. An ac voltage with an 
amplitude of 0.05 V and a frequency of 33 Hz was ap- 
plied to the strain gauges. Compared with our previous 
data taken by a dc method, the signal-to-noise ratio was 
markedly improved. 12 

When external pressure was applied, the pressure cell 
(a NiCrAl - BeCu hybrid piston cylinder cell) was heated 
to about 320 K to reduce inhomogeneity arising from 
the solidification of Daphne oil (see discussion below for 
details). The pressure at low temperatures was deter- 
mined by the superconducting transition temperature of 
indium. 

2.2 P — T phase diagram of Daphne oil 7373 
Inhomogeneity will exert a harmful effect upon SmS. 

Since the pressure inhomogeneity is induced by the solid- 
ification of the pressure-transmitting medium, it is im- 
portant to know the "frozen-melted" phase diagram of 
Daphne oil 7373. We start with the construction of its 
P — T phase diagram using SmS as a test material. 

Figure 1(a) shows the temperature dependence of the 
thermal expansion AL/L(T) of SmS at 10 and 13.5 kbar. 
The anomaly observed at approximately 250 K is due to 
the freezing of Daphne oil. This enables us to determine 
the "solidification pressure" or equivalently the "melting 
temperature" of Daphne oil. 

The solidification pressure is shown in Fig. 1(c) as a 
function of temperature. When Daphne oil is kept at 
room temperature, it does not solidify even at 20 kbar. 
Reflecting this nature of Daphne oil, the electrical re- 
sistivity of a strain gauge, which was attached to the 
reference sample Cu, shows no anomaly up to about 
22 kbar, as shown in Fig. 1(b). By contrast, the same 
measurement using Fluorinert as a pressure-transmitting 
medium shows an anomaly at approximately 10 kbar, 
which we ascribe to the solidification of Fluorinert. 13 

The P — T phase diagram of Daphne oil tells us that 
pressure should be applied at the highest temperature 
possible. This is the reason why we heated the pressure 
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Fig. 1. (Color online) (a) Temperature dependence of thermal ex- 
pansion AL/L(T) of SmS at 10.0 and 13.5 kbar measured during 
heating. The anomalies denoted by arrows are due to the solidi- 
fication of Daphne oil 7373. (b) Pressure dependence of the rela- 
tive electrical resistivity of a strain gauge at room temperature. 
Here, the gauge was attached to the reference sample Cu, and 
Fluorinerts FC70 and FC77 and Daphne oil 7373 were used as 
pressure-transmitting media. The arrows show anomalies due to 
the solidification of the pressure-transmitting media, (c) Freez- 
ing temperature of Daphne oil 7373. The circles and triangles 
were deduced from the temperature and pressure dependences 
of thermal expansion and electrical resistivity shown in (a) and 
(b), respectively. The broken line is a visual guide. 



cell when applying pressure. This is also effective in the 
case of Daphne oil 7474. 14 

2.3 P — T phase diagram of golden SmS 

Figure 2(a) shows the temperature dependence of the 
thermal expansion coefficient cx(T) obtained from the 
thermal expansion AL/L(T) given in the inset via the 
relation a(T) = (dL/dT)/L. At ambient pressure (black 
phase), the sample length L, proportional to the volume 
V for the cubic crystal, decreases monotonically with de- 
creasing temperature (see the inset). At 3.1 kbar (golden 
phase), AL/L(T) exhibits a strikingly different behavior 
from that of the black phase; the volume increases with 
decreasing temperature below about 130 K. Correspond- 
ing to this unusual behavior in the AL/L curves of the 
golden phase, a(T) exhibits a broad minimum (~ —4.5 x 
10~ 5 ) at a temperature To. Note that this minimum value 
is almost independent of pressure. The peak structure 
is described by a Schottky model and the characteristic 
temperature T decreases with pressure. These results 
are consistent with our previous reports, 4, 12 indicating 
that there is no distinction between the data obtained for 
Daphne oil and Fluorinert in the pressure range tested 
(P < 16.2 kbar), although the pressure range exceeds the 
solidification pressure of Fluorinert (^10 kbar). 

At pressures above P C 2 ~ 19 kbar, a kink anomaly 
appears in the dilatation AL/L(T) curve at a low tem- 
perature. Corresponding to that, a sharp anomaly occurs 
in the a(T) curve, which signals a Neel transition, T^. 
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Fig. 2. (Color online) (a) Temperature dependences of thermal 
expansion coefficient a(T) of SmS at selected pressures up to 
21.2 kbar. To and Tjv are defined as the temperatures at which a 
shallow minimum and a positive, rather sharp peak appear, re- 
spectively. The inset shows the temperature dependence of AL/L 
at the same pressure as that of o(T). Data at different pressures 
are shifted for clarity, (b) Temperature dependences of thermal 
expansion coefficient at selected pressures between 16.2 and 19.1 
kbar. The inset shows o(T) along the broken line in the main 
frame as a function of pressure. The region sandwiched by the 
two broken lines denotes a two-phase mixture region. 

The positive sign of the anomaly implies that Tjv in- 
creases with pressure, consistent with the present results 
and previous data of specific heat and magnetic suscep- 
tibility. 5 ' 8 

Figure 2(b) shows the temperature dependences of 
thermal expansion coefficient at different pressures be- 
tween 16.2 and 19.1 kbar. Interestingly, we observe a 
peak with a positive sign; the peak is superposed on the 
broad minimum situated at about 10 K and grows with 
increasing pressure along the broken line. It is reason- 
able to ascribe the observation to a two-phase mixture 
consisting of the broad minimum and positive peak aris- 
ing from the paramagnetic and antiferromagnetic states, 
respectively. 

In the inset of Fig. 2(b), we plot a along the broken 
line. This is interpreted as the volume fraction of the an- 
tiferromagnetic phase, which steeply increases at a pres- 
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Fig. 3. (Color online) Temperature dependences of the compress- 
ibility k^(T) at selected pressures. The broad peak below 17.8 
kbar and the sharp peak above 19.6 kbar correspond to To and 
Tjv deduced from a(T), respectively. Data at different pressures 
are shifted for clarity. See text for details. 

sure interval AP ~ 1.8 kbar between 17.3 and 19.1 kbar. 

Figure 3 shows the temperature dependences of 
the compressibility Kl(T) at selected pressures, which 
were obtained using the relation kl — ~{AL/T|p 2 — 
AT/T|p 1 }/(P 2 — -Pi)- At low pressures, we only observe 
a broad peak corresponding to the anomaly of a(T) at 
To. As pressure increases, the peak structure becomes 
prominent. At pressures higher than P C 2, e.g., 19.6 kbar, 
a sharp peak appears at the Neel temperature. 

From these anomalies, we construct the P — T phase 
diagram of golden SmS, as shown in Fig. 4. The bro- 
ken line through the open and closed circles denotes a 
crossover dividing the paramagnetic phase into the low- 
and high-temperature regions. The solid line through the 
open and closed squares indicates the second-order Neel 
transition, and the double solid line through the triangles 
corresponds to the first-order Neel transition, as shown 
below. The large circle denotes the TCP that separates 
the first-order transition from the second-order transi- 
tion. The hatched region indicates the two-phase mix- 
ture region mentioned above. Interestingly, note that the 
crossover line Tq(P) (deduced from a(T j) appears to ter- 
minate at the TCP, although the data points near TCP 
deduced from Kl deviate from it presumably owing to 
pressure inhomogeneity. According to Matsubayashi et 
al., Tq is the characteristic temperature of the formation 
of a pseudogap. 4 

To estimate dilatation as a function of pressure 
AL/L(P), instead of temperature, we need to know the 
absolute value of the compressibility k (or equivalently 
the bulk modulus B) at a certain temperature lower than 
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Fig. 4. (Color online) P — T phase diagram of golden SmS. Circles 
and squares were deduced from a(T) and k^(T) shown in Figs. 
2 and 3, respectively. The broken line denotes a crossover at To, 
and the solid and double solid lines correspond to the second- 
order and first-order Neel transitions, respectively. The hatched 
region at approximately P c 2 shows the two-phase mixture region 
that is possibly caused by sample inhomogeneity. The inset shows 
the pressure dependences of the thermal expansion AL/L(P) at 
5 and 12 K. Note that the AL/L(P) curve shows a jump at 
approximately 18.5 kbar. This anomaly defines the first-order 
Neel transition plotted in the phase diagram. The large circle 
shows the tricritical point (Pt, Tt). 




Fig. 5. (Color online) Pressure dependence of electrical resistivity 
at 2 K. It shows a maximum at approximately 15.5 kbar, which 
does not coincide with P c 2 (= 18.1 ± 0.8 kbar). Note that the 
data points above P c 2 are correspond to the antiferromagnetic 
phase. 



the volume change is ascribed to the two-phase mixture 
arising from the inhomogeneity of pressure and/or the 
sample. 

According to our thermal expansion study of CeRhlns 
and UGe2, the pressure gradient of Daphne oil 7373 is 
estimated at less than 0.2 kbar. 17 ' 18 This is much smaller 
than the pressure interval AP mentioned above. We con- 
sider that the wide transient region is attributed to sam- 
ple inhomogeneity rather than to pressure inhomogene- 
ity. As mentioned in Introduction, the pressure-induced 
black-to-golden phase transition at P c \ is accompanied 
by a large volume change that amounts to about 10%. 
Therefore, it is likely that some inhomogeneity is induced 
in the sample during the black-to-golden phase transi- 
tion. 



the melting temperature of Daphne oil, e.g., at 150 K (see 
Fig. lc). However, we only have available data taken at 
room temperature. Therefore, we assume that B(T — 
150 K) = B(T — 300 K). 15 Any other choice of B does 
not change the following conclusion. The thus-evaluated 
AL/L{P) at a low temperature is shown in the inset of 
Fig. 4. At 5 K, AL/L(P) shows a steplike, albeit ob- 
scured by some inhomogeneity (see below), anomaly at 
approximately 18 kbar. By contrast, no jump is visible at 
12 K; we only observe the P-linear behavior of AL/L(P). 
This means that the TCP lies between 5 and 12 K. From 
a detailed study of the T dependence of the jump AV, 16 
we determined the position of the TCP on the P — T 
plane to be (P t , T t ) = (19.0 ± 0.8 kbar, 10.5 ± 0.5 K), 
as shown in Fig. 4. Note that the compressibility exhibits 
a very large anomaly near the TCP (Fig. 3). We also de- 
termine P C 2 by the extrapolation of the first-order line 
to zero temperature as P C 2 — 18.1 ± 0.8 kbar. 

Strictly speaking, the volume change observed is not a 
jump but a rapid variation at a narrow pressure interval. 
Note that this transient region between 17.5 and 19.1 
kbar almost coincides with the AP mentioned above (see 
the inset of Fig. 2(b)). Therefore, the broadened nature of 



2.4 Pressure dependence of residual resistivity 

Let us discuss in more detail the volume jump in con- 
junction with valence change. In general, a Sm 2+ ion has 
a greater ionic radius than a Sm 3+ ion. Furthermore, a 
conduction electron, which can be produced as a result 
of the promotion of a 4f electron from a Sm 2+ ion into a 
conduction band, tends to expand its wave function so as 
to obtain a kinetic energy gain, resulting in a reduction 
in crystal volume. As a result, it is expected that a semi- 
conducting phase composed of Sm 2+ ions has a larger 
volume than a metallic phase composd of Sm 3+ ions. In 
the present case, therefore, it is very likely that the av- 
erage volume is larger in the paramagnetic phase with 
the mixed valence than in the antiferromagnetic phase 
with the metallic conductivity that mainly consists of 
Sm 3+ ions. This suggests that valence change occurs at 
the border of antiferromagnetism. 

According to theory, a valence change gives rise to an 
enhancement of residual resistivity. 19 This is because im- 
purity potential is enhanced by a many-body effect in the 
vicinity of a quantum critical point of valence transition 
at which the correlation length of valence fluctuations di- 
verges. Therefore, we expect that the residual resistivity 
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of golden SmS shows a maximum at approximately P c2 . 

To check this, we measured electrical resistivity un- 
der pressure. The detailed results will be published else- 
where. In Fig. 5, we show only the pressure dependence of 
the electrical resistivity taken at 2 K, which is tentatively 
regarded as the residual resistivity po- It is clearly seen 
that po initially increases with pressure and then passes 
through a maximum before steeply decreasing, consistent 
with previous reports. 20 ' 21 It is possible to speculate that 
the initial increase in p is ascribed to the decrease in mo- 
bility from the increase in carrier concentration at 300 K 
with pressure. 22 

As is expected, po shows a maximum near P c2 . How- 
ever, we note that the maximum pressure does not ex- 
actly coincide with P c2 . There are two possible explana- 
tions for this deviation. First, the phase transition at P c i 
is of the first-order; therefore, the theoretical model men- 
tioned above may not be applicable to the present case. 
Second, the existence of the transient region gives rise to 
not only the broadening of the peak structure but also 
the shift in peak position. It is unclear at present which 
is more probable. Thus, we need further investigation. 

3. Conclusions 

We first constructed a phase diagram of Daphne oil 
7373. The freezing temperature at 25 kbar was estimated 
as about 300 K. Therefore, we pressurized the cell at 320 
K to avoid the inhomogeneity induced by the solidifica- 
tion of the pressure-transmitting medium. 

We measured the thermal expansion of golden SmS at 
pressures of up to 22.7 kbar, and observed a broad peak 
at the characteristic temperature T and a sharp peak 
at the Neel temperature Tm- From these anomalies, we 
constructed the P-T phase diagram. It consists of three 
characteristic lines: the crossover line corresponding to 
To, and the second-order and first-order phase lines cor- 
responding to T/v. These lines appear to emerge from the 
tricritical point, located at (T t , P t ) ~ (19.0 ± 0.8 kbar, 
10.5 ± 0.5 K) on the pressure vs temperature plane. A 
linear extrapolation of the first-order phase line provides 
the critical pressure P C 2 = 18.1 ± 0.8 kbar. In the vicin- 
ity of P c2 , there is a two-phase- mixture region, in which 
both the paramagnetic and antiferromagnetic phases are 
mixed owing to sample inhomogeneity 

We observed a volume jump at P C 2 where the quantum 
phase transition between the paramagnetic and antiferro- 
magnetic phases occurs. We measured residual resistivity 
as a function of pressure, and found that it shows a max- 
imum at a pressure near P C 2 . Considering the correlation 
of volume with valence, we pointed out the possible va- 
lence jump at the border of antiferromagnetism. 
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